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Abstract 

The nature of the phase transition for the XY stacked triangular antifer- 
romagnet (STA) is a controversial subject at present. The field theoretical 
renormalization group (RG) in three dimensions predicts a first order tran- 
sition. This prediction disagrees with Monte Carlo (MC) simulations which 
favor a new universality class or a tricritical transition. We simulate by the 
Monte Carlo method two models derived from the STA by imposing the con- 
straint of local rigidity which should have the same critical behavior as the 
original model. A strong first order transition is found. Following Zumbach 
we analyze the second order transition observed in MC studies as due to a 
fixed point in the complex plane. We review the experimental results in order 
to clarify the different critical behavior observed. 

P.A.C.S. numbers:05.70.Fh, 64.60.Cn, 75.10.-b 

I. INTRODUCTION 

Phase transitions of frustrated spin systems have been extensively studied during the 
last decade (for reviews see In particular the nature of the phase transition of the 

stacked triangular antiferromagnet (STA) with XY spins interacting via nearest-neighbor 
bonds has been extensively investigated At high temperatures the symmetry group of 

this system is 0(2) Z2 whereas at low temperatures this symmetry is completely broken. 
The Ising symmetry Z2 has its origin in the non coUinearity of the spins in the ground state 
which can be classified as chirality plus or minus. For non frustrated systems the symmetry 
group in the high temperature region is simply 0(2) and this difference in symmetry between 
frustrated and non frustrated spin systems should lead to a different critical behavior. Bailin 
and Garel 0, using the renormalization group (RG), proved that there is no stable fixed 
point close to space dimension d = 4 and they concluded that the transition is of first order. 
Extending the RG technique to a iV component spin system, that is using 4 — e expansion 
to first order in e (two- loops) and expanding also in 1/N, Kawamura suggested a new 
universality class linked to the chirality for the transition with the exponents given by Monte 
Carlo simulations [0]. With the same technique in d = 4 — e and in three dimensions to 
three-loops, it was shown later that the transition for N = 2 must be of first order P,^. 
Further Monte Carlo studies have confirmed that the exponents for STA-system (see 



table 0) are different from the ones of the standard 0{N) universahty class (given in table 
Vl). Simulations seem to favor the concept of a new chiral universality class or tricritical 
behavior. However, Plumer and Mailhot |T0| used different exchange constants for spins 



along the c-axis and spins in the triangular planes, so that the hexagonal STA-system is 
quasi one-dimensional. They concluded that the transition is weakly first order. 

There are two principal groups of magnetic materials which can be modeled by our 
system. The first group are Hexagonal perovskites ABX3, which are quasi one dimensional 
systems which however order at low temperatures and have a planar anisotropy so that the 
spins are in the XF-plane. The most studied examples are CsMnBrs pTj^T|, RbMnBrs 
CsVBrg g, CsVClg |§, and CsCuCls For a review see ref. The results 



of the first four compounds are compatible with second order transitions with exponents 
more or less in agreement with the MC simulations: for example u = 0.50(1) in MC and 
between 0.54(3) and 0.57(3) experimentally (for details see Tables 1). However, the specific 
heat measurement of CsCuCls indicates a cross over to first order in zero magnetic field . 

The second group are helimagnetic systems. Since the angle between the spins can 
be different from 120° for the STA-structure without changing the critical behavior ||11|| , 
helimagnetic rare earths (see ref. for Ho [^9| -|56[|, for Dy ||57|-|63| and for Tb 0-^]) could 
also be analyzed by the STA-model. For helimagnets the critical behavior is quite varied 
(see the review for Ho and Dy ref. [0). Essentially three types of results exist: one in favor 
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of a O (4) class ||5^ , |55|j6l| , |62[| , another in favor of a new universality class 
and a third class which favor first order transition [p0| , p7| , ^ . See Tables for details. 



With these results a definite answer cannot be given about the order of transition. In section 
V we will come back to this point. 

In order to check the results of the renormalization group studies we have studied 
the STAR and the Stiefel model |T^. The first is derived from the STA model by imposing 
the constraint that in each triangle the sum of the spins is zero at all temperatures. The 
modes removed are irrelevant for the RG and the two models STA and STAR should be 
in the same universality class, provided such a class exists. As is explained in the next 
section the Stiefel model we use for the simulation is connected to the STAR model. Each 
cell of three spins plus constraint is equivalent to a system of dihedral, i.e. an ensemble 
of two perpendicular vectors. Neighboring pairs of vectors interact ferromagnetically, but 
only vectors of the same kind. Since these two systems have the same number of degrees of 
freedom, they should belong to the same universality class. We think that the two models 
are closer to the RG studies than the original stacked triangle antiferromagnet. We can 
therefore check the predictions of the RG and gain an understanding of the difficulties one 
has with the results of Monte Carlo simulations and measurements in the critical region. 

In section II we present the two models. The simulations and the details of the finite 
size scaling analysis for a first order transition are explained in section HI. The results are 
shown in section IV. Discussion and conclusions are in section V. 



II. MODEL AND SIMULATION 
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A. The STAR model 



Starting from the stacked triangular antiferromagnet (STA) we take the simplest Hamil- 
tonian with one antiferromagnetic interaction constant J > 

H = JJ2S..S,, (1) 

where are two component classical vectors of unit length. The sum runs over all nearest 
neighbor pairs, that is the spin Sj has six nearest neighbor spins in the same XF-plane and 
two in Z-direction in adjacent planes. The ground state is characterized by a planar spin 
configuration with three spins on each triangle forming a 120° structure with either positive 
or negative chirality (see Fig. 1). The ground state degeneracy is thus twofold like the Ising 
symmetry, in addition to the continuous degeneracy due to global rotations. 

In the RG theory one uses the concept of local rigidity which means that the sum of 
three spins Si, S2, S3 on the corners of a triangle is set to zero 

Si + S2 + S3 = . (2) 

In this theory the local fluctuations violating this constraint become modes with a gap. Thus 
they do not contribute to the critical behavior and can be neglected |TB|. The constraint in 



(0) used for an ordinary collinear antiferromagnet with two spins instead of three eliminates 
also one degree of freedom so that only one is left which means that the critical behavior of 
antiferromagnet is the same as that of a ferromagnet. In our case we are left with two degrees 
of freedom. We choose one spin direction and then have one more choice for the chirality, 
that is the direction of the second spin could be chosen clockwise or counter-clockwise with 
respect to the first spin. 

In order to impose local rigidity for the MC simulation, we first partition the lattice 
into interacting triangles which do not have common corners. This can be done as follows. 
In each XF-plane one selects in a row one " supertriangle" . Then one finds two nearest 
supertriangles which do not share a common corner in the row (they are separated by two 
head-up and three head-down triangles). This process is repeated. Then one takes the next 
row until all rows in the XF-planes are filled with triangles which do not share corners, see 
Fig. 1. All spins are then located on the supertriangles and each spin belongs to only one 
supertriangle. Local rigidity means that the three spins in each supertriangle form a 120° 
structure. Only in the ground state all the supertriangle have the same orientation. At finite 
temperature local rigidity means that there are no local fluctuations within a supertriangle, 
but fluctuations between supertriangles are allowed. 

The MC updating procedure for the state of the supertriangles is made as follows. At a 
supertriangle, we take a new random orientation for one of its three spins; next we choose 
a second spin so as to form a ±120° angle with the first spin. For the orientation of the 
third no freedom is left. The interaction energy between the spins of this supertriangle with 
the spin of the neighboring supertriangles is calculated in the usual way and we follow the 
standard Metropolis algorithm to update one supertriangle after the other. 

We consider L * L * Lz systems , where L * L is the size of the planes, and = 2L/3 
the number of planes. L must be a multiple of three so that no frustration occurs because 
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of periodic boundary conditions in the XF-planes. Simulations have been done for systems 
sizes with L = 12, 18, 24, 30, 36. 

The order parameter M used in the calculation is 

M = jrj:\Ms\, (3) 

s=l 

where Mg {s = 1, 2, 3) is the s-th sublattice magnetization and = L'^L^ is the total number 
of the lattice sites. This definition corresponds to the one for the ordinary antiferromagnet 
with only two sublattices. Instead of alternating signs for the coUinear case the non collinear 
staggered magnetization is obtained by making a rotation of +120*^ (—120°) for the second 
(third) magnetization before summing over the three sublattice magnetizations. 
The chirality k is defined in the usual way 

^ [S] X + X + Sf X Si] , (4) 



3^3 
1 



where the summation is over all supertriangles and A^' = A^/3 is their number. The chirality 
Ki of one triangle is parallel to the Z-axis and equal to ±1. 



B. The Stiefel model 



The Stiefel model can be derived from the STAR model |T^. We give the main points. 
In each elementary cell an orthonormal basis 

ea(^); a = 1,2 (6) 

is defined, where i is the superlattice index. Each spin located in the cell can be represented 
in this basis 

S. = E^a(Oea(0 • (7) 
a 

If we put this expression into the Hamiltonian (|l]) we obtain a new Hamiltonian with inter- 
actions between the orthogonal vectors ea{x): 

H = JY.[ ei(0 ■ ei(j) + e2(z) ■ e^U) ] . (8) 

The interaction can be chosen negative (or ferromagnetic) and the sum is for simplicity 
over nearest neighbor pairs of a simple cubic lattice instead of a hexagonal lattice since the 
new spins (see Fig. 2 and 3) are no longer frustrated. The chirality k, for the Stiefel 
model is defined as 

1 

« = ^lEiei(i) X e2(i)| . (9) 
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The Hamiltonian (Jg) is similar to the one of the Ashkin- Teller model [jT4[. Indeed we 
can give this Hamiltonian a form which is close to it. The interaction energy of two nearest 
neighbors (ij) with opposite chirality (|^) is zero and with the same chirality it is 2ei(i)ei(j) 
(see Fig. 3). Therefore the Hamiltonian can be written as 

H = JJ2{'^ + ^^^j)^^^j (10) 

where a = ±1 is an Ising spin representing the chirality and Sj is an XF-spin. In Hamilto- 
nian of Ashkin and Teller only Ising spins appear. 

Despite the fact that the Stiefel model is extensively studied, especially in two dimensions, 
no clear picture emerged. The problem is to know whether there are two transitions, an 



Ising and a Kosterlitz-Thouless transition, or only one transition of a new type [^,^. In 
three dimensions it has been shown that there is only one transition Here the 

problem is to determine the order of the transition. 

The procedure of MC procedure is as follows. At each site one takes a new random 
orientation for the first vector and chooses for the second vector a perpendicular direction 
(we have two choices: ±90°, the Ising degrees of freedom). We have two degrees of freedom 
the same number as for the STAR model. The interaction energy between this dihedral and 
its neighbors is calculated. If it is lower than the energy of the old state, the new state 
is accepted. Otherwise, it is accepted only with a probability according to the standard 
Metropolis algorithm. It is possible to use a cluster MC algorithm |T^, but in the case of a 



strong first order transition there is no reduction of the critical slowing down |T7|. Periodic 
boundary conditions are used. 

Systems with L = 12, 15, 18, 21, 24 have been simulated. To compare with the size L 
of the STA or the STAR model, we must multiply L by a/3. One supertriangle or triangle 
contains three spins and is represented by one site in the Stiefel model. So we obtain 
equivalent sizes from 20 to 40. 

The order parameter M for this model is 

1 2 

M = y\MJ, (11) 

2iV ' ' ' ^ ^ 

where Mg (s = 1, 2) is the magnetization for the vectors over all sites and = is the 
total number of sites. 



C. Definitions, liistogram metliods and finite-size scaling 

We use in this work the histogram MC technique developed by Ferrenberg and Swendsen 
T^,[TP|. The histogram for the energy Pt{E) is very useful for identifing a first order tran- 



sition. Also the data obtained by simulation at To can be used to obtain thermodynamic 
quantities at temperature T close to Tq. Since the energy spectrum of a Heisenberg spin 
system is continuous, the data list obtained from a simulation is basically a histogram with 
one entry per energy value. In order to use the histogram method efficiently, we divided the 
energy range < by 10 000 bins. We have verified that we obtain the same results, with 
our precision, for 30 000 bins. 
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The critical slowing down in a first order transition is greater than in a second order 
transition because of energy barriers, and thus the time of the simulation, to go from one 
state to another grows exponentially with the size of the lattice. For this reason we restricted 
our simulations to systems not too large to have good enough statistics. In each simulation, 
at least 2 millions (3 millions for the greater sizes) measurements were made after enough 
Metropolis updating (500 000) were carried out to reach equilibration. 

For each temperature T we calculate the following quantities 

C - «E'>-<E>^) 

N(<M^>-<M>^) 
^ = ' 

X. - ^^<^'>-<^>') , (14) 

where M is the order parameter, C the specific heat per site, x the magnetic susceptibility 
per site, V the fourth order energy cumulant, < ... > means the thermal average. 

The finite size scaling (FSS) for a first order transition has been extensively studied 
72l- A first order transition should be identified by the following properties: 



a) Pt{E) has a double peak. 

b) The maximum of the specific heat C and the susceptibilities x ^^^d Xk are proportional 
to the volume L'^. 

c) The minimum of the fourth order energy cumulant V varies as: 

V = V* + bL-'^, (16) 

where V* is different from 2/3. 

d) The temperatures T{L) at which the quantities C, x or Xk have a maximum should 
vary as: 

T{L) = T, + aL-'^. (17) 
All this conditions will be verified for our systems. 



III. RESULTS 

In Fig. 4 the specific heat C of the STAR model is plotted as function of the temperature 
for various sizes (we note that the maximum is 30 times as large as the usual value of STA 
which is a sign of a first order transition). The value of the maximum as function of the 
volume is shown in Fig. 5 for C, x aiid Xk- We note that the maxima vary like the volume 
except for the smaller sizes where further corrections are important. 
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In Fig. 6 the same quantities are shown for the Stiefel modeL In all cases the maxima 
vary for the greater sizes proportional to the volume as they should for a first order transition. 

In Fig. 7 and 8 we have plotted V as function of T for different sizes L for the STAR 
and the Stiefel model respectively. We can see that V does not tend to 2/3 (for a second 
order transition) but to a value V* < 2/3. This value is calculated by fiting the minimum 
with ([T6|) . As result we yield for the STAR model 

V* = 0.652(2) (18) 

and for the Stiefel model 

V* = 0.625(3) . (19) 

The Fig. 9 and 10 show the energy distribution for different sizes at different tempera- 
tures for the STAR and Stiefel model. The double peaks observed, even for a very small L, 
indicate a strong first order transition. With increasing sizes, these two peaks are separated 
by a region of zero probability, indicating a discontinuity of the energy at the transition. 
We estimate the correlation length by 1/3 of the first size where the two peaks are well 
separated by a region of zero probability. This is our estimate of the distance needed for 
two phases to coexist. This method yields the correct answer in the case of Potts models. 
We obtain ~ 12 for the STAR model and ~ 9 for the Stiefel model (see Fig. 9 and 10). 

To obtain the critical temperature we can use (0). The results are 

Tc = 2.2990(5) (20) 

for the STAR model and 



Te = 2.4428(4) 



(21) 



and for the Stiefel model. We have sizable corrections for the small systems. Comparing the 
last results with those of Kunz and Zumbach |T2[ for the case V2,2 similar to ours we agree 
with their result Tc=2.445. 

Our results show clearly the first order transition for the STAR and Stiefel model. So we 
confirm the indication given in for the Stiefel model where in high temperature region 
u was determined not too far away from 1/3 which is the value for a first order transition. 



IV. DISCUSSION 

We have shown that the STAR and the Stiefel model have first order transitions. These 
models are equivalent to the STA in the RG theory, because the constraint of local rigidity 
is not relevant in the transition region. Consequently the phase transition for the original 
triangular antiferromagnet STA must also be of first order, and this result holds generally 
for all systems with a breakdown of symmetry of type 0(2) ® Z2. How can we reconcile our 
results with those of the MC of the STA model and the experimental results? 

For the MC simulation of the original frustrated spin system the second order transition 
is an effect of the finite system size according to Zumbach's analysis of "almost second order 



phase transitions" |25|. The main point of this analysis is that the stable fixed point F^, 
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known to exist only for the number of components > N^, moves into the unphysical 
complex plane when N < N^. In our case N = 2, the estimation for iVc is iVc = 3.91 |§] 
and second order could never occur. Nevertheless this complex fixed point has a large basin 
of attraction and mimics a behavior of a real fixed point. Only if the system is very large, 
i.e. if L > where is the largest correlation length, the transition will appear of first 
order. The phenomenon of a crossover between a second order to a first order transition 
is not so uncommon. An extreme case is the Potts model in two dimensions with q = 5 
components, where the transition is known to be of first order [p^. The MC gives always 



a second order transition with critical exponents of an instable fixed point due to the 
enormous correlation length ,^o- 

If the Hamiltonian is a sum of two terms, one interaction is of the Heisenberg symmetry 
and the other favors Ising symmetry: H = Hneisenberg + Hjsing and if Hjsing < Hneisenberg, 
one has a crossover between a region of Heisenberg type to Ising behavior close to T^. If the 
system size is too small we will only see the region controlled by the Heisenberg fixed point 
and therefore obtain the exponents of the Heisenberg magnet. In a sense we have the same 
situation if we replace the Heisenberg fixed point by the Zumbach fixed point. However, 
there is an essential difference: for the fixed point in the complex plane, one has to modify 
the scaling relation 



7/z/ = 2 — ?7 + c 



(22) 



by a constant c different from zero p5|. If the fixed point is real c must be zero. 

We can use this relation as a criterion for real or complex fixed points. In three dimensions 
rj is usually small, that is ~ 0.03 (see table [V^) , but it must be positive [^] and therefore 
7/z/ < 2 for a genuine second order phase transition. If the ratio 7/z/ > 2 the fixed point 
must be complex. The correction c in the scaling law (^) will depend on the distance of the 
fixed point from real space. Therefore one expects that c will be greater for the XY {N = 2) 
case than for the Heisenberg (A^ = 3) case. For frustrated Heisenberg systems it will be 
difficult to find out whether rjefj = 77 — c is negative. 

With the relation 7/z/ = 2 — 1] we obtain rj = —0.06 using the results of MC simulation 



and 77 = —0.16 from experimental values of Ho ||5^. Thus the fixed point is in the complex 
plan and the second order transition observed in the XY systems has only an "almost second 
order" character. 

The effect of imposing local rigidity obviously forces the system to stay away from the 
region of influence of the complex fixed point F^. and thus permits to "see" the true first 
order behavior. Introducing larger coupling constants for inter-plane interactions than for 
intra-plane interactions |jTO| seems to have a similar effect as the local rigidity constraint. 



A further remark concerning MC simulations and the fixed point Fc in the complex plane; 



Zumbach has shown that the finite size scaling in this case does not hold. The FSS 
ansatz for the free energy should be replaced by 



f = L~''g[L/^,C2HL)] 



(23) 



where 5^ is a function of L/^, but also of InL. The constant C2 is proportional to the constant 
c of (P^. This constant is small and therefore the correction to FSS. Indeed if we take for 
the true value of r] the value of the ferromagnetic case r] ~ 0.03 (see table [V| ) the value of 
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c will be at most equal to c ~ 0.03 + 0.06 ~ 0.1 (see discussion above) and if we compare 
with 2 — 1] in (|2^) this gives an error of 5 %. However, small but not negligible corrections 
to the standard FSS could explain the differences in MC simulations obtained with different 
methods (see table 0) and also some of the differences in the experimental values (tables 

We will discuss now the experiments in the light of the concepts used. In order to see 
the first order region the temperature resolution is the limiting factor not the finite size. 
However, they are linked through to oc ^q^^'^ with u ~ 0.5 found by MC and depending 
on the materials studied. The temperature distance to = {T — Tc)/Tc could be too small to 
be measurable. 

The experiments on CsMnBra RbMnBrg and CsVBra Q (D give expo- 



nents compatible with those of MC on STA and a second order transition (see table | and 
We can interpret this result by the fact that the systems are under the influence of a 
complex fixed point and to is too small to observe a first order transition. 
The case CsCuCls (table D is different since the authors observe a crossover from a 
second order region with exponents compatible with MC on STA for 10"^ < t < 5.10"^ to 
a region of first order transition for 5.10"^ < t < 5.10"^. For t < to ~ 10"'^ one seems to 
observe the true first order region. 

Helimagnetic rare earth metals are more complicated as already discussed in the intro- 
duction (see also tables ^^V]). The results compatible with those of the MC on STA for 
Ho [g-g (table HI), Dy |3|,|9],|g (table |T|) and Tb |4|-||] (table ^ can be interpreted 
as before: the systems are under the influence of Fc. The first order transition for Ho ||5^ 



and Dy p7| , |58| is due to the fact that the measurements were done in the first order region 
near the critical temperature. The values of the exponent (3 ~ 0.39 in the case of Ho and 
Dy (see table pil-[m|) are not compatible with those found by MC {P ~ 0.25). This fact 
can be explained by the presence of a second length scale in the critical fluctuations near 
Tc related to random strain fields which are localized at or near the sample surface ||52[| . 
Thus the critical exponent f3 measured is of this second length. However the result of j3 for 
Tb (table ITVl ) shows values compatible with MC but it has been proved that this second 
length is present also in Tb Further measurements to determine (3 should help in the 
interpretation. 

We have tried to give a general picture of the phase transition of frustrated XY spins 
where the breakdown of symmetry is of type 0(2) ®Z2. We have shown that this transition is 
first order but usually not seen because of the presence of a fixed point in the complex plane. 
One way to observe that the behavior is really driven by such a fixed point is the existence 
of a negative value for t] in the Monte Carlo simulations and experimental systems. The 
method used here is to impose local rigidity. This constraint does not change the behavior 
of the system but permits the system to stay outside the region of influence of the complex 
fixed point. So we can rely on the renormalization group study and the true behavior is first 
order. There is no "new chiral universality class" in strict sense for our system (A^ = 2). 
Another possibility discussed in the literature is that the transition is influenced by the 
presence of topological defaults which are not visible in the continuum formulation of the 
RG (for the presence of topological defaults in Stiefel model see [0). 

From our experience with the frustrated XF-model we conclude that the true first order 
transition for the frustrated Heisenberg model cannot be reached in MC simulations. The 
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experimental situation should be similar. However, due to presence of the anisotropies we 
never reach the Heisenberg first order region but will have a crossover to the Ising or XY 
region . 
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TABLES 



Crystal 



method 



ref. 



CsMnBrs 



Neutron 



0.22(2) 



CsMnBr'^ 



Neutron 



0.25(1) 



CsMnBr^ 



Neutron 



0.24(2) 



CsMnBrs 



Neutron 



0.21(2) 



1.01(8) 



CsMnBrs 



Neutron 



1.10(5) 



CsMnBrs 



Calorimetry 



0.39(9) 



CsMnBrs 



Calorimetry 



0.40(5) 



RbMnBrs 



Neutron 



0.28(2) 



RbMnBrs 



Calorimetry 



0.22-0.42 



CsCuCl. 



Neutron 



0.25(2) 



CsCuCls 



Calorimetry 



0.35(5) if 10~^ < t < 5.10- 



CsCuCls 



Calorimetry 



> 0.6 if 5.10-5 < t < 5.10" 



TABLE I. Experimental values of critical exponents for compound AXB3 



Crystal 
Thermal expansion 



ref. 



a 



1 order 



Calorimetry 



0.27(2) 



Calorimetry 



0.10-0.22 



Neutron 



1.14(10) 



Neutron 



0.3(1) 



1.24(15) 



Neutron 



0.39(2) 



Neutron 



0.39(4) 



X-ray 



0.37(1) 



X-ray 



0.39(4) 



TABLE II. Experimental values of critical exponents for Holmium 
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J.1 1 

method 


ref. 


a 


(3 


7 


V 


inermal expansion 






1 order 




Calorimetry 






1^* order 




Caiorimetry 




r nT 


n 1 o / o \ 

U.18(8j 












6^ 










Neutron 




53|| 






1.05(7) 


0.57(5 


Neutron 




Fl 




0.38(2) 






Neutron 




61| 




0.38(3) 






Neutron 




62| 




0.39(1) 






Mossbauer 




63| 




0.335(10) 







TABLE III. Experimental values of critical exponents for dysprosium 



method 


ref. 


a 


/3 


7 


V 


Calorimetry 


11 


0.20(3) 








Neutron 


H 




0.25(1) 






Neutron 






0.23(4) 






Neutron 


n 








0.5 


X-ray 






0.21(2) 







TABLE IV. Experimental values of critical exponents for Terbium 



ref. 



a 



0.34(6) 



0.253(10) 



1.13(5) 



0.54(2) 



-0.09(8) 



0.45(2) 



Ik 



0.77(5) 



0.55(2) 



0.46(10) 



0.24(2) 



1.03(4) 



0.50(1) 



-0.06(4) 



0.38(2) 



0.90(9) 



0.55(1) 



i 



0.43(10) 



0.48(2) 



TABLE V. Critical exponents by Monte Carlo for 0(2) ® Z2. ""^calculated by 7/z/ = 2 — ry. The 
first result [0] comes from a study at high and low temperature and uses of FSS. The second 0] 



uses the Binder parameter |24| to find and uses the FSS, the third Q uses the maxima in FSS 
region. The results /3k, 7k, Vk are the exponents for the chirality. 
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symmetry 
Z2 
0(2) 
0(3) 
0(4) 
0(5) 
0(6) 



a 

0.107 ' 
-0.010 ' 
-0.117' 
-0.213' 
-0.297' 
-0.370 ' 



/3 
0.327 
0.348 ' 
0.366 ' 
0.382 ' 
0.396 ' 
0.407' 



1_ 

1.239 
1.315 
1.386 
1.449 
1.506 
1.556 



0.631 
0.670 
0.706 
0.738 
0.766 
0.790 



TABLE VI. Critical exponents for the ferromagnetic systems calculated by RG |26| 
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FIGURES 




FIG. 1. Ground state configuration for tlie 
STA and the STAR model. The chirahty of 
each triangle is indicated by + or — . The other 
ground state configuration with opposite chiral- 
ity is obtained by reversing all spins. The super- 
triangles are dark. 



a) 




b) 

FIG. 3. a. Dihedral and their interaction. 
ei{i) interacts only with ei{j) not with 62(7) 
which interacts with 62 (^). b. Two dihedral with 
opposite chirality. The energy of the interaction 
is equal to cos(7r — 9) + cos{9) = 0. 



L\ 



L\ 



FIG. 2. The dihedral are drawn at the cen- 
ter of each elementary supertriangle. 
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FIG. 4. Specific heat for the STAR model 
for various sizes as function of temperature. 



FIG. 6. Maxima of the specific heat C, the 
susceptibility x (magnetization) and Xk (chiral- 
ity) in function of the volume V = for the 
Stiefel model. 
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0.655 



»^ =0.652(2) 
2.30 2.31 



2.32 



FIG. 5. Maxima of the specific heat C, the 
susceptibility x (magnetization) and Xk (chiral- 
ity) in function of the volume V = !?■ ■ for 
the STAR model. 



FIG. 7. V as function of temperature for 
different sizes for the STAR model. The . . . line 
indicate the value of F* = 2/3 for the second 
order transition. The broken line is our estimate 
of V* at the critical temperature Tc = 2.2990. 
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0.67 



0.66 



0.65 



0.64 



0.63 



0.62 
2.439 





1 


Z.=24 ^/^'''^^ 




1 


' / / 

\ / / L=2^ 1 

/ 






\ )\^\ / / 
\ V V\ / 


/ L=15 






— ^Z.=12 


l^*=0.625(3) 





2.449 



2.459 



FIG. 8. V as function of the tempera- 
ture for different sizes for the Stiefel model. 
The . . . Hne indicate the value of V^* = 2 /3 for 
the second order transition. The broken line is 
our estimate of V* at the critical temperature 
Tr = 2.4428. 



FIG. 9. Energy histogram P{E) as function 
of energy per site E for the STAR model for var- 
ious sizes L at different temperatures of simula- 
tion: Ti2 = 2.3136, Tig = 2.3065, T24 = 2.3020, 
Tso = 2.30085, Tge = 2.29968. 
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FIG. 10. Energy histogram P{E) as func- 
tion of energy per site E for the Stiefel model 

for various sizes L at different temperatures 
of simulation: Tu = 2.4495, T15 = 2.4468, 
Tis = 2.44525, T21 = 2.44425, T24 = 2.44377. 
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